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Abstract 

Langmuir turbulence is an archetype of wave turbulence in plasma physics. By means of ID- IV 
Vlasov-Poisson simulations, we show that coherent structures, called Langmuir cavitons, are generated 
by the long time evolution of Langmuir weak turbulence, thus illustrating the breakdown of a weak 
turbulence regime. These structures correspond to an equilibrium between the pressure forces and 
the ponderomotive force resulting from high frequency Langmuir oscillations. Langmuir cavitons are 
typical features of strong Langmuir turbulence expected to be generated at high energy and to saturate 
when Langmuir energy is of the order of the plasma thermal energy. Despite this wide-spread belief, 
here we observe that cavitons, emerging from weak Langmuir turbulence evolution, saturate at much 
lower energies. We show that these Langmuir coherent structures are characterized by a much larger 
length scale with respect to the Debye length. This gives evidence that "large" and "shallow" stable 
cavitons should be seen in space plasma observations. The transition toward strong turbulence is 
shown to be a consequence of an initial weak turbulent inverse cascade. Finally, the effective equation 
of state for ion acoustic oscillations is tested numerically from the kinetic model. 



1 Introduction 

The nonlinear evolution of waves is usually classi- 
fied in terms of "weak" and "strong" turbulence. 
The meaning of these terms is not well established 
although the major conceptual difference between 
weak and strong turbulence is the presence of a 
characteristic dimensionless parameter e charac- 
terizing the level of nonlinearity, as the ratio be- 
tween a typical linear (or dispersive) time (the elec- 
tron plasma period in Langmuir turbulence) and 
the nonlinear time, depending on the wave ampli- 
tudes. Turbulence is then considered to be "weak" 
if e ^ 1; otherwise it is "strong". A second as- 
sumption is usually made: weak turbulence results 
from a superposition of finite, but weak amplitude 
waves, obeying the linear dispersion relation, but 
with randomly distributed phases. Starting with 
(almost) random phases, such randomness must be 
preserved over the nonlinear evolution time. Two 
key points are decisive for the validity of Random 
Phase Approximation (RPA): (a) the amplitude of 



the fields and (b) the bandwidth of the phenom- 
ena under consideration [1] [51 [3]. At finite but 
small energy, the nonlinear dynamics is described 
by three-waves or four-waves interactions in the 
RPA through the kinetic wave equations, whereas 
at higher energies intermittency dominates the dy- 
namics, through the apparition of coherent struc- 
tures. 

To summarize, weak turbulence is based on the 
following points: (i) linear dispersion still holds; 
(ii) statistical homogeneity in space holds; (iii) 
wave particle interaction is described by resonant 
quasilinear theory. On the other hand, strong 
turbulence is based on the following: (i) linear 
dispersion no longer holds; (ii) strong statistical 
inhomogeneity in space; (iii) wave particle interac- 
tion become complex because of wave (in density 
cavities) and particle trapping. 

Only a few attempts have been made to compare 
weak turbulence theory with numerical results. 



The three dimensional equations for capillary 
water waves has been solved numerically [1] and 
a Zakharov & Filonenko power-law spectrum [5] 
was observed. Weak turbulence theory may be 
valid in some spectral range but fails in others 
[HI- Recently, the problem of breakdown of weak 
turbulence by intermittent events associated 
with coherent structures has been addressed 
for different models [71 [6] and an illustration 
of the modulational instability influence on the 
breakdown of weak turbulence resulting from an 
inverse cascade has been given for the nonlinear 
Schrdinger equation [5]. However, these results 
have been established for fluid models, i.e. without 
kinetic effects, whose role in the phase randomiz- 
ing is unknown. Therefore, the validity of weak 
turbulence theory in a full kinetic regime is still 
an open problem to be investigated by numerical 
simulations. 



1.1 Langmuir turbulence 

We hereafter concentrate on the specific case of 
electrostatic Langmuir turbulence, an archetype of 
wave turbulence in plasma physics. In the weak 
turbulence regime, the 3-wave (resp. 4-wave) evo- 
lution dominates for large (resp. small) wave num- 
bers fcL > ki\ii (resp. fci < kMi) through the decay 
(resp. modulational) instability, which typically 
transfers the L-wave energy towards smaller (resp. 
larger) wave vectors. The transition wavenum- 
ber is knji^o =1/3 Cs/vth,e with Aj^ the Debye 
length, Cs the ion sound speed and Vth,e the elec- 
tron thermal velocity [5]. Weak Langmuir turbu- 
lence is mainly driven by the electrostatic decay 
of a Langmuir wave (hereafter L-wave) into an- 
other L-wave and an ion acoustic wave (hereafter 
lA), while several processes associated to the pon- 
deromotive force, as modulational instability, os- 
cillating two stream instability, soliton formation, 
Langmuir collapse, are considered to dominate the 
strong turbulence dynamics [HE]. In the specific 
case of Langmuir turbulence, the turbulence is usu- 
ally considered "strong" for an electric-to- thermal 
energy ratio larger than the electron-to-ion mass 
ratio: W = cqE^ /nkBT > rrie/mi. 

Langmuir cavitons are localized electric fields os- 
cillating at the plasma frequency self-consistently 
associated to density cavities OD] with a large range 



of variations, 0.01 < Sn/n < 0.80 [TT]. These struc- 
tures have been observed in the laboratory in sev- 
eral plasma experiments |12j while in natural space 
plasmas only in active ionospheric experiments |13j . 
Cavitons at "high" energies have also been inten- 
sively studied through numerical experiments of the 
Zakharov equations [TH [TSl dni E] and the Vlasov- 
Poisson equations [TSl [HI [20] . For moderate forc- 
ing, it has been shown that weak turbulence and 
strong turbulence features can coexist [H] . 

The Langmuir cavitons canonical scenario, 
starting from an initial long wavelength Langmuir 
spectrum, can be summarized as follows: (a) 
parametric cascading leading to Langmuir conden- 
sation and increase in wave intensity; (b) cavity 
formation {e^E'^ /2nkBT > k'^X'jy); (c) caviton 
collapse {k — )■ oo); (d) particles acceleration and 
cavity emission of ion sound waves; (e) start of a 
" caviton nucleation cycle" . 

By considering here the limit of ID Langmuir 
turbulence, an archetype of wave turbulence in 
plasma physics, we show the transition from weak 
to strong turbulence through the formation of co- 
herent structures (cavitons), independently of the 
initial level of coherence of Langmuir oscillations. 
The resulting cavitons may saturate at "low" en- 
ergy levels (electric energy orders of magnitude 
lower than thermal energy) generally considered 
to belong to the weak Langmuir turbulence regime 
and remain then relatively stable. We obtained a 
power law governing the relation between the scale 
length of the structures and their energy to be di- 
rectly tested on space plasma data. Finally, the 
nonlinear time scale needed to reach the strong tur- 
bulence regime is found to scale as the inverse of the 
initial Langmuir energy. 

These results are obtained with a kinetic descrip- 
tion of ID electrostatic plasma. To our knowledge, 
no similar work on the breakdown of weak turbu- 
lence has been done using a fluid model (Zakharov 
equations) . 

2 Model 

We solve the IDIV Vlasov-Poisson system of equa- 
tions for the electron and proton distribution func- 
tion /e , fp and the self-consistent electric potential 
and electric field, and E. All equations are nor- 
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malized by using electron quantities, the electron 
charge, mass and thermal velocity, e, rrie and Vth.e, 
the plasma frequency ujpe, the Debye length and 
a characteristic density and electric field, fie and 

E = ■meVth,e^pe/e. 

We define W = 0.5 x as the electric energy 
density normalized to the electron kinetic energy. 
Then, the dimensionless Vlasov equations read: 



dfe 

dt 



^dfe 

dx 



- (E + E.^t)'^ ^ 
ov 



ot ox fi ou 



(1) 



(2) 



where v and u are the electron and ion velocity, 
/i — m^/mp — 1/1836 the electron-to-proton mass 
ratio and E^xt an "external" driver acting on the 
electrons only that can be switched on or off dur- 
ing the runs. The details of the external forcing 
can be found in Appendix 1 of Ref. [21 . Finally, 
the Vlasov equations are self-consistently coupled 
to the Poisson equation. 



= / fedv - / fpdu ; E = - 



dx 



(3) 



We use a numerical box of length = 5000 A„ and 
a velocity range —5 < v/vth,e < +5 for electrons 
and —5 < u/uthA < +5 protons. The spatial and 
velocity mesh grid is dx = X^, dv = 0.04 Vth.e and 
du = 0.04 Uth,i where Uth,i is the proton thermal 
velocity. Periodic boundary conditions are used 
in the spatial direction. The initial electron and 
proton velocity distributions are Maxwellian with 
equal temperatures, Tp = T^, a typical condition in 
solar wind plasma, so that lA fluctuations are effi- 
ciently damped out. We add at t = a density ran- 
dom noise in the wavelength range 30 < A < 500. 

We consider two different initial conditions, (i) 
coherent and (ii) incoherent L- waves with an initial 
electric energy level varying in the range 10""* < 
W < 10^^, corresponding to the transition from 
weak to strong turbulence regime. In case (i), the 
external driver Ef,xt, oscillating at ujpe, excites a 
monochromatic L-wave with wavelength A^^ and 
propagating in one direction only. The forcing is 
switched off when the L-wave reaches the desired 
amplitude i?^. The initial wave then evolves self 
consistently according to the Vlasov-Poisson sys- 
tem. The corresponding Langmuir energy ranges in 



the interval 2 • 10~* < W < 10~^, here again cor- 
responding to the transition from weak to strong 
Langmuir turbulence regime. Finally, the initial 
random density noise is 5n/n = 10^^. In case (ii), 
we exclude the forcing from the beginning and let 
the system to evolve starting with an electron den- 
sity random noise corresponding to a flat spectrum 
in the electric field. The r.m.s. electron density 
amplitude is 10~^ < Sue/rie < 10~^, correspond- 
ing to a r.m.s. electric field 10^^ ^ E < 10^^, i.e. 
an energy range 2 • 10"'* < W < 2 ■ 10"^. 

We use a kinetic model to fully take into account 
all processes, as the wave-particle interactions, 
that could limit the development of Langmuir 
turbulence by extracting electric energy and con- 
verting it into kinetic energy. The phase velocity 
V(j, of the initial L- waves is much larger than Vth,e, 
so that the Landau damping of Langmuir waves is 
inefficient. On the other hand, Tp = Tg (a common 
situation in space plasmas like the solar wind) 
so that the generated lA waves are efficiently 
Landau damped during the transient part of the 
simulation, corresponding to a weak turbulence 
evolution. This kinetic damping could limit, first 
the Langmuir cascade, second, the generation 
of density fluctuations that are a seed for the 
subsequent generation of density inhomogeneities 
observed in the asymptotic part of the simulation, 
corresponding to the strong turbulence evolution. 



3 Numerical results 

We flrst discuss the results of a simulation from 
the first set of numerical experiments, starting 
with a monochromatic L-wave with E^^ = 0.06 
and Aj^ = 100 A^ . The time evolution of the 
electric energy and the ion density are shown in 
Fig. [l] The L-wave first undergoes electrostatic 
parametric instability (Langmuir electrostatic 
decay) during the period 3 • 10** < t < 6 ■ W^. 
In the mean time, lA fiuctuations are generated. 
Then, ion cavities start to form at t ~ 10^. These 
cavities are filled by electric energy in the form 
of an electrostatic field oscillating at the plasma 
frequency. The ion density fluctuations 6n/n and 
the Langmuir energy Wj^ are shown in Fig. [2j 
top panel, (blue and black lines respectively) for 
t > 1.5 X 10^. The ion density fluctuations as well 
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Figure 1: Long time evolution of electric energy 
(top panel) and ion density Sn/n (bottom panel) 
in the {x,t) plane, starting from a monochromatic 
Langmuir wave. Note in the bottom panel the gen- 
eration of (i) ion acoustic waves from the Langmuir 
electrostatic decay between 3 x 10"* < t < 6 x 10"*; 
(ii) ion cavities filled with electric energy for t > 
10^. 



Figure 2: Top panel: ion cavitons with the as- 
sociated ion density fluctuations Sn/n (blue line) 
and electric energy (black line). Bottom panel: 
The pressure and the ponderomotive force (red and 
black line, respectively). The ion density fluc- 
tuations as well as the envelop of the Langmuir 
electric energy remain practically constant from 
t > 1.5 X 10^ 



as the envelope of the Langmuir electric energy 
remain practically constant from t > 1.5 x 10^ to 
the end of the numerical simulation. The cavities 
result from an equilibrium between the total pres- 
sure, — V(Pe + Pi) and the ponderomotive force, 
—e"^ / {iniiUJp^) dxE"^, associated to high frequency 
Langmuir oscillations (see Fig. [2] bottom panel, red 
and black line, respectively). These coherent struc- 
tures, identified as "cavitons", arc the signature of 
a transition to a strong Langmuir turbulent regime. 



We have performed many simulations starting 
with (i) coherent and (ii) incoherent initial L-waves. 
In all cases, the system freely evolves until it even- 
tually jumps to a strong turbulence state charac- 
terized by the presence of Langmuir cavitons. The 
formation of cavitons does not need a high level of 
L-waves if the system evolves "long enough". In 
case (i), as discussed, parametric decay occurs first 
and saturates in the first part of the simulation. 



Then cavitons are formed and remain stable un- 
til the end of the simulation. In case (ii) the ini- 
tial electron density fluctuation self-organises into 
a large spectrum of Langmuir noise. The "sea" 
of L-waves that fills the simulation box then col- 
lapses into stable low energy cavitons. In all cases 
the nonlinear structures are similar, indicating that 
the asymptotic behavior of Langmuir turbulence 
is independent from the initial level of coherence 
of Langmuir oscillations. As expected, regardless 
of the initialization of the turbulence, the larger 
the initial L-waves, the sooner Langmuir cavitons 
are generated. For each caviton, we have measured 
the depth of the ion density hole 6n/n, the length 
of the structure L (width of the ion cavity at the 
height of 1 /e of the maximum depth) and the max- 
imum Langmuir electric energy Wj^ that sustains 
the density hole. The results are shown in Fig. [3] 
on a Wj^ interval ranging over three decades. It 
is worth noticing that Langmuir cavitons are ob- 
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served also at remarkably low electric energy val- 
ues, ~ 10"'^. The depth of the ion cavity 
5n/n is of the order of the Langmuir electric en- 
ergy density (as expected for high energy cavitons) 
over three decades: 

6n/n = {Sn/n)„ (4) 

with {Sn/n)^ = 0.28 ±0.06 and a = 1.13 ±0.06 the 
fitting parameters. Fit results are given with their 
respective 3cr errors. The width of the cavitons also 
scales on the Langmuir electric density energy: 

L = L„ W^f (5) 

with L„ = 18 ± 4 and /3 = -0.47 ± 0.05. As 
expected, Langmuir cavitons have a larger scale 
length when the Langmuir electric energy is lower. 
A new unexpected and important result is the gen- 
eration of stable Langmuir caviton with a scale 
length of many hundred of Debye lengths. Finally, 
the power law between the depth and the length of 
the ion cavities associated with the cavitons results 
as: 

L = L„ {5n/ny (6) 
with = 10 ± 3 and 7 = -0.42 ± 0.05. 
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Figure 3: Depth expressed in relative density Sn/n 
(top panel) and width L expressed in Debye length 
(bottom panel) of cavitons measured in the sim- 
ulations according to the associated Langmuir en- 
ergy . Each diamond represents a single cavi- 
ton. The line shows the power law fit. 



4 Discussion 

We hereafter discuss, first, the difference between 
the cavitons observed in our simulations and 
the Langmuir solitons also associated to strong 
Langmuir turbulence. Then we investigate the 
transition from weak to strong turbulence and its 
observed timescale. 



4.1 Langmuir solitons and cavitons 

Interestingly, the scaling laws obtained here for 
the Langmuir cavitons, Eq. |4]|6] are similar to 
those of Langmuir solitons [5^]. In particular, 
the typical depth of the ion cavities is of the 
order of the Langmuir electric energy, obtained 
in the simulation over three orders of magnitude. 
We recall that Langmuir cavitons, or Langmuir 
standing solitons, are propagating structures main- 
taining their shape through the balance between 
dispersive and nonlinear effects, while cavitons are 



standing structures with electric oscillations that 
self-consistently sustain the density cavity of which 
they are eigenmodes. 



4.2 Transition from weak to strong 
turbulence 

In order to identify the mechanism responsible for 
the transition from weak to strong turbulence, we 
have studied the temporal evolution of the electric 
and ion density spectra. In the first part of the ex- 
periments, successive Langmuir decay instabilities 
generate smaller and smaller wave numbers (3- wave 
inverse cascade characteristic of weak turbulence). 

The 3-wave cascade spontaneously develops for 
the first series of numerical experiments starting 
from a monochromatic L-wave. This is shown in 
Fig. |4] (initial wave vector — 27r/50) where we 
draw the development of the Langmuir turbulence 
(top panel) and the excitation of daughter S-waves 
(bottom panel) driven by the successive decays of 
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Figure 4: Evolution of the electric field and 
ion density spectrum (top and bottom pan- 
els respectively). Black dashed lines: ex- 
pected wavenumber from 3- wave cascade in- 
teractions. Blue dashed line: transition 
wavenumber kj^^j between decay and mod- 
ulational instabilities. 

the L-waves. The L-wave cascade foUows the path 
in Fourier space expected for such 3-wave cascade 
process (the theoretical wave numbers are shown 
in dotted lines). As previously pointed out [5T], 
the decay product S-waves are generated over a 
finite band of wave numbers and do not survive 
outside the coupling region with L-waves (bottom 
panel) because of strong Landau damping. Note 
that there is no cascade on the S-wave: instead 
the harmonics of the S-waves are temporarily ex- 
cited as long as a pump L-wave injects energy into 
a fundamental S-wave, but the high kinetic damp- 
ing prevents any further development of an inertial 
range on these S-waves. 

Interestingly, the development of the L-wave cas- 
cade allows to identify isothermal electrons and adi- 
abatic ions as the effective equation of state for ion 
acoustic fluctuations. Here the ion sound speed Cs 
is given by = {jeTe + jiTi)/mi, where 7e = 51 
and 7i = 5/3. A slightly different definition of Cs 
would lead to an incorrect path of the cascade in 



Fourier space. This effective equation of state, well 
known in the linear regime with Te >> T^, is thus 
extended to ~ Ti in the weak turbulence regime. 
This result is confirmed by a scattered plot of the 
ion pressure vs ion density during the weak turbu- 
lent stage of the simulation (not shown here). 

Once this weak turbulent cascade has devel- 
oped, cavitons are systematically generated (ap- 
parition of a large spectrum in both panels) when 
fc^Aj-, ~ 10^^ corresponding to a critical wave vec- 
tor kMi^o = 1/3 Cs/vth^e (blue dashed line in 
Fig. |4| for which the growth rate of the modula- 
tional instability, a known precursor of strong tur- 
bulence, overcomes the growth rate of the decay 
instability [9 . This is confirmed by other simula- 
tions starting with different initial wave vectors. 

To sum up, the weak turbulent inverse cas- 
cade brings the fluctuations to smaller wave vec- 
tors, where different nonlinear processes generate 
the coherent structures typical of a strong turbu- 
lent regime. In this picture, the strong turbulence 
regime is the natural asymptotic limit of weak tur- 
bulence. 

Finally, a series of experiments has been carried 
out with an initial wave vector fcz, < kMi, in order 
to directly trigger the modulational instability in- 
stead of the decay instability. No L-wave cascade 
is observed, but similar cavitons are indeed gen- 
erated, in agreement with the previously exposed 
picture. 

The detailed mechanism for the generation of 
the cavitons is still unclear, even if it is known that 
cavitons can be by-products of the development 
of the modulational instability. This point will be 
investigated in future works. 



4.3 Evaluation of the nonlinear 
timescale 

To complete this previous picture, we show in Fig. [5] 
the characteristic nonlinear time scale t^^^ corre- 
sponding to the formation time of the first caviton 
in a simulation. We see that r„j^ scales as the in- 
verse of the initial Langmuir energy W,^ -^^^ : 

^ r„,o Wl^^^ (7) 

with r„^o = 178 ± 4 and 77 = -1.0 ± 0.1. Since 
the growth rate for Langmuir decay scales as 
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Figure 5: The nonlinear time corresponding to the 
formation of Langmuir cavitons for different initial 
Langmuir energies W^^ .^^.^ . Diamond: simulation 
results, blue line: fit. 

W^i,,i„it [12 ! ^NL is interpreted as the sum of 
a succession of Langmuir decay time scales to 
reach the critical wave vector kc for which the 
modulational instability becomes dominant and 
can generate the observed coherent structures 
(cavitons). Therefore, in a plasma like the solar 
wind, since in the L- waves are typically observed 
in the range 10~^ < W°'"' < 10~^, the time scale 
to reach a strong turbulence regime assuming 
Wf'' - 10^2 ^ouid bg abo^t iQ (^ith fpe 10 
kHz). 



5 Conclusions 

In summary, we have shown that L-wave turbu- 
lence is expected to breakdown for times larger 
than a typical nonlinear time scale, in agreement 
with Refs. [IllS]. The distinction between weak and 
strong turbulence thus looses part of its significa- 
tion. These results are independent from the ini- 
tial level of coherence. The formation of Langmuir 
cavitons appears to be the asymptotic limit of weak 
turbulence, as a result of an inverse Langmuir cas- 
cade. The evolution of the electric spectrum shows 
that the weak turbulence regime brings itself the 
disturbances to large enough wavelengths for which 
the condensation of coherent structures becomes ef- 
ficient, thus enabling the transition to strong tur- 
bulent features. The kinetic description allows to 
identify the effective equation of state for ion acous- 
tic fluctuations in the weak turbulence regime of an 
electrostatic plasma with equal electron and proton 



temperatures. It is given by isothermal electrons 
and adiabatic ions; this extends the result known 
in the limit Te >> T^. Furthermore, electrostatic 
coherent structures of typical width much greater 
that a few Debye lengths are generated by the long 
time evolution of an initial relatively moderate am- 
plitude turbulence. These results have been ob- 
tained for a ID kinetic description of electrostatic 
plasmas and should be extended to multidimen- 
sional studies in order to test the stability of the 
observed coherent structures. 

The breakdown of weak Langmuir turbulence 
and the existence of large coherent structures can 
have an important impact on the interpretation 
of space plasma data. The authors are thus 
confident that these new insights in Langmuir 
turbulence may encourage the space physics 
community to revisit the admitted conclusion 
that strong turbulent Langmuir structures are 
formed at too high energies to be relevant in space 
plasma environments. This last result should be 
directly tested on waveforms data in space plasma 
environments. 
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